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The dilute A3 model is a solvable IRF (interaction-round-a-face) model with three 
local states and adjacency conditions encoded by the Dynkin diagram of the Lie 
algebra A3. It can be regarded as a solvable version of a critical Ising model in a 
magnetic field. One therefore expects the scaling limit to be governed by Zamolod- 
chikov's integrable perturbation of the c = 1/2 conformal field theory. We perform 
a detailed numerical investigation of the solutions of the Bethe ansatz equation for 
the off-critical model. Our results agree perfectly with the predicted values for the 
lowest masses of the stable particles and support the assumptions on the nature of 
the Bethe ansatz solutions which enter crucially in a recent thermodynamic Bethe 
ansatz calculation of the factorized scattering matrix. 



1 Introduction 

Despite the fact that the two-dimensional Ising model in a magnetic field has 
defied an analytic solution to this day, relatively much is known about the 
properties of this model. In particular, a celebrated result of Zamolodchikov^''^ 
shows that a perturbation of the c = 1/2 conformal field theory^''* (which corre- 
sponds to the critical Ising model) with the spin density operator of conformal 
dimensions (A, A) = (1/16, 1/16) leads to an integrable quantum field theory 
which contains eight massive particles with factorized scattering, the par- 
ticle masses and the ^'-matrix elements being related to the exceptional Lie 
algebra Eg. 

In an appropriate scaling limit, the Ising model in a magnetic field should 
be described by this integrable field theory, and the predictions for the lowest 
mass ratios have been examined numerically by several authors. As these 
results rely on relatively small size transfer matrix calculations or on Monte 
Carlo simulations, they provide rather crude checks for the lowest mass ratios 
only. 

The recent construction ^ solvable version of the Ising model in a 

field, the so-called dilute A3 model, gives us the opportunity to have a closer 
look at this scaling limit. By a thermodynamic Bethe ansatz calculation, 
it has indeed been possible to recover Zamolodchikov's result for the mass ra- 
tios and the ^'-matrix. This calculation relies on assumptions concerning the 
string structure of the solutions of the Bethe ansatz equations for the dilute 
A3 model. However, a first numerical check shows that the largest eigenvalues 
of the transfer matrix at criticality do not conform with the predicted string 
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structures. This observation has motivated us to perform an extensive numer- 
ical investigation of the corresponding off-critical Bethe ansatz equations, to 
obtain precise estimates of some of the mass ratios on the one hand and in 
order to check the assumptions on the structure of its solutions on the other. 

This article is organized as follows. We commence by introducing the 
dilute A3 model. In Sec. 3, the Bethe ansatz equations for the dilute A3 model 
are given. Subsequently, we discuss the critical behaviour of the model, and in 
Sec. 5 the scaling limit is defined. In Sec. 6, we present our numerical results 
and we conclude in Sec. 7. 

2 The dilute A3 model 

The so-called dilute A-D-E models are a series of solvable lattice models (in 
the sense of commuting transfer matrices) which have been constructed re- 
cently. 16,13 xhey are critical IRF (interaction-round-a-face) models with local 
states and adjacency conditions encoded by Dynkin diagrams of the simply 
laced Lie algebras. The dilute A^ models can be extended to off-critical so- 
lutions of the Yang-Baxter equation with Boltzmann weights parametrized in 
terms of elliptic functions of the spectral parameter. These models are of par- 
ticular interest for an odd number of local states L. In contrast to the critical 
case, the off-critical weights of these models do not respect the Z2-symmetry 
of the underlying Dynkin diagram. This means that the solution describes 
the critical system in the presence of a symmetry-breaking field rather than a 
temperature-like deviation from criticality. 

Here, we are only interested in the dilute A3 model, which is a three-state 
model built on the effective adjacency diagram depicted in Fig. 1. In an allowed 
configuration on the square lattice, all pairs of states on neighbouring lattice 
sites have to be directly connected on the adjacency diagram. In other words, 
the states '-|-' and ' — ' cannot occur on adjacent lattice sites in any allowed 
configuration. The model is defined by associating Boltzmann weights 

d c 



a b 

to each configuration of states a,b,c,d around the elementary plaquettes of the 
square lattice. As the explicit form of the weights is rather cumbersome, we 
refrain from quoting these here and follow the notation of Refs. 15, 17 and 18. 
Let us only mention that besides their dependence on the spectral parameter 
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Figure 1: Effective adjacency diagram of the dilute A3 model. 



u, the weights depend on the elliptic nome q (with q = corresponding to the 
critical case). For fixed q, they satisfy the Yang-Baxter equation 



e d e d 




ah ah 



which implies ^® that the row transfer matrices T[u) form a one-parameter 
family of commuting matrices, i.e., [T(m), T(t;)] = 0. 



3 Bethe Ansatz Solution 

The dilute models can be solved by a Bethe ansatz. eigenvalues 
A(m) of the row transfer matrix T[u) (with periodic boundary conditions) are 
given by 

^ ( ^i(m-2A) ^i(m-3A) \^ -jX- ^i(M-Mj+A) 
~ V ^i(2A)^i(3A) ) ^i(M-Mi-A) 



??i(2A) ??i(3A) J ^}-^{)i(u-Uj-\){)i(u-Uj-2\) 



1 ( di{u)di{u-\)\^ ^ di{u-Uj-A\) 



where N denotes the number of faces in the row, and 



A = 57r/16, w = exp(i7r£/4), < m < 3A , (4) 
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Figure 2: Sketch of typical arrangements of Bethe ansatz roots Vj in the complex plane. 



where £ G {1,2,3}. The Bethe ansatz roots Uj are determined by the Bethe 
ansatz equations 

f di{uj-\) \^ - _ TT ^i("j-"fc-2A) di{uj-Uk+\) , , 

with i = 1,2,...,7V. 

Defining Vj = iuj, the largest eigenvalue Aq(m) (in the sector I = V) 
corresponds to a set of real Bethe ansatz roots Vj. In general, the solutions 
will however be complex, but for large system size N they typically arrange 
into so-called strings (subsets of roots with approximately the same real part) 
as shown schematically in Fig. 2, where the horizontal line represents the real 
axis. It should be noted that the 'string type' of a particular solution might 
well depend on the variable q, an example of this behaviour will be given below. 



4 Critical Behaviour 



As already mentioned, the dilute A3 model is critical for g = 0. In the regime 
under consideration, see Eq. 4, the critical point is of Ising type and the model 
is therefore described by a c = 1/2 conformal field theory with scalar primary 
fields of conformal dimensions (A, A) with A = A = (corresponding to the 
identity operator), A = A = 1/16 (the spin density operator) and A = A = 
1/2 (the energy density operator). Each of these primary fields gives rise to an 
infinite tower of descendant fields with conformal dimensions (A + r, A + r), 
where r and r are positive integers, which together form the product of the 
corresponding irreducible representations of the Virasoro algebra. 

As a consequence of conformal invariance, the appropriately scaled spectral 
gaps 

= £ log(Ar/Ar) (6) 

in the critical limit (g = 0, — ?> c») have the form A + r + A + r with 
the aforementioned values. Here, we denoted by K^p the eigenvalues of the 
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Table 1: Finite-size approximants for central charge c and smallest scaling dimensions Xj. 



N 


c 


Xl 


X2 


xs 


X4 


X5 


X6 


X7 


10 


0.499681 


0.1250802 


1.0145 


2.196 


3.148 


4.272 


4.417 


5.449 


20 


0.499920 


0.1250200 


1.0035 


2.142 


3.034 


4.062 


4.191 


5.098 


50 


0.499987 


0.1250032 


1.0006 


2.128 


3.005 


4.010 


4.135 


5.015 


75 


0.499994 


0.1250014 


1.0003 


2.126 


3.002 


4.004 


4.130 


5.007 


100 


0.499997 


0.1250008 


1.0001 


2.126 


3.001 


4.002 


4.128 


5.004 


oo 


1/2 


1/8 


1 


2 1/8 


3 


4 


4 1/8 
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transfer matrix T(3A/2) at the isotropic point u = 3A/2 for q = Q, Aq being 
the largest eigenvalue. Moreover, the central charge c = 1/2 manifests itself in 
the finite-size corrections of the largest eigenvalue 

-log(A(°)) = Nfo + ^ +o{N-'). (7) 

In Table 1, numerical approximants for the smallest scaling dimensions Xj of 
zero momentum operators (i.e., A + r = A + r) are presented for several system 
sizes N. These are obtained from Eq. 6, where the eigenvalues Aj°' have been 
calculated by solving the Bethe ansatz equations numerically, see Sec. 6 for 
details. Furthermore, Table 1 also includes finite-size approximants for the 
central charge c = 1/2 obtained by neglecting the o{N~^) correction terms in 
Eq. 7, where the known analytic expression ^® for the specific free energy /o 
has been used. 

5 Scaling Limit 

The additional parameter q (which has the physical interpretation of a mag- 
netic field) allows us to approach the critical point in different ways, resulting 
in a scaling limit characterized by the scaling variable jj, 

^ = qN^'l\ (8) 

i.e., the limits g — and N ^ oo are taken simultaneously keeping jj, constant. 
Here, = corresponds to the critical limit discussed above, jj, = oo io the 
massive limit where Zamolodchikov's Eg field theory applies, see Fig. 3. The 
appropriately scaled spectral gaps 

/, = log(Ao/A^.) (9) 

in the scaling limit become functions of the scaling variable jj, alone. 
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Figure 3: Schematic view of the scahng hmit in the variables q and A^. 



In the massive limit, the ratios rj = supposedly approach the par- 

ticle mass ratios of the corresponding massive field theory. The mass ratios of 
the eight stable particles are determined by the entries of the Perron- Frobenius 
eigenvector of the Cartan matrix of the Lie algebra Eg and (ordered by mag- 
nitude) are given by 



nii/m 


= 1 








1712 /m 


= 2 


cos 




= 1.618034... 


nis/m 


= 2 


COS 


V30) 


= 1.989044... 


1714 /m 


= 4 


COS 


7r/5) cos(77r/30) 


= 2.404867... 


1715 /m 


= 4 


COS 


7r/5) cos(27r/15) 


= 2.956295... 


nie/m 


= 4 


COS 


7r/5) cos(7r/30) 


= 3.218340... 


nir/m 


= 8 


COS' 


(7r/5) cos(77r/30) 


= 3.891157... 


nis/m 


= 8 


COS' 


(7r/5) cos(27r/15) 


= 4.783386... 



(10) 



where the corresponding enumeration of nodes in the Dynkin diagram of Eg is 
shown in Fig. 4. 




Figure 4: Dynkin diagram of Eg with nodes labelled according to the masses of Eq. 10. 
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6 Numerical Results 

For the numerical treatment of the Bethe ansatz equations we employed a mod- 
ified Newton method We took the following course of action: As a first step, 
the eigenvalues of the transfer matrix were obtained by direct diagonalization 
for small system sizes (up to = 10). Secondly, the corresponding solutions 
of the Bethe ansatz equations were identified and generalized to larger N. 
Finally, these solutions were then followed as the variable q is varied. 

Since we are mainly interested in the masses of single particle states, we 
restricted ourselves to the zero momentum sector, i.e. to translation invariant 
eigenvectors of the transfer matrix. The numerical data of Table 1 confirm that 
we correctly identified the states with scaling dimensions a; < 5 in the critical 
limit. Apart from the largest eigenvalue, which corresponds to a set of A^ real 
Bethe ansatz roots Vj for all values of q (0 < q < 1), the root arrangements of 
all the other eigenvalues change several times as q is varied. A simple case is 
shown in Fig. 6, where the Bethe ansatz roots for the second largest eigenvalue 
(sector £ = 2) are given for N = 2. Comparing different system sizes, one 
finds that these changes actually occur at approximatively constant values of 
the scaling parameter so they do enter in the scaling limit. Of course, the 
corresponding finite-size eigenvalues are analytic functions of q always. 

As mentioned above, the thermodynamic Bethe ansatz calculation ^® relies 
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Figure 6: Ratios of scaled gaps as functions of the scaling parameter ^. The corresponding 
scaling dimensions (A + r, A + r) are indicated. 



on assumptions on the string structure of the Bethe ansatz solutions. For 
the eigenvalues under consideration, we find that for sufficiently large scaling 
variable /j, these assumptions are indeed correct, except for one eigenvalue 
(corresponding toA = A=l/2 and r = r = I (x = 3), respectively to the 
mass 1114) where one might have to go to much larger values of the scaling 
variable than we did. Therefore, the apparent confiict mentioned in Sec. 1, 
between the string types which are observed at criticality and those which 
enter the thermodynamic Bethe ansatz approach, appears to be solved. 

In Fig. 6, the dependence of the ratios rj = of the scaled gaps fj 

(Eq. 9) on the scaling variable /j, is shown. In this plot, we used data from 
system sizes N = 50, = 75 and N = 100, they clearly scale as expected. This 
is also supported by the numerical data for three values of scaling parameters 
/J, given in Table 2. The convergence to the limit values is quite obvious, 
and one recognizes four non-trivial mass ratios as well as two ratios belonging 
to two-particle states (with masses 2mi and nii + 1112, respectively). Also, 
among the scaling functions considered we observe one intersection (between 
the functions labelled by conformal dimensions (3/2, 3/2) (mass 1114) and (2, 2) 



8 



Table 2: Ratios of scaled gaps for three values of the scaling parameter fi. 



A* 


N 


ri 


r2 




r4 


r5 


re 




10 


1.618749 


1.983943 


2.403558 


2.083461 


2.616910 


2.910107 




20 


1.618037 


1.983923 


2.403649 


2.081017 


2.617335 


2.930942 


10 


50 


1.618024 


1.983977 


2.403797 


2.079804 


2.617468 


2.934096 




75 


1.618024 


1.983984 


2.403815 


2.079550 


2.617482 


2.934583 




100 


1.618025 


1.983990 


2.403834 


2.079551 


2.617497 


2.934617 




10 


1.623692 


1.989446 


2.412441 


2.017536 


2.623691 


2.958322 




20 


1.618144 


1.987688 


2.405009 


2.025316 


2.618144 


2.954218 


20 


50 


1.618034 


1.987597 


2.404866 


2.026212 


2.618034 


2.954289 




75 


1.618034 


1.987591 


2.404865 


2.026286 


2.618034 


2.954306 




100 


1.618034 


1.987589 


2.404865 


2.026315 


2.618034 


2.954311 




10 


1.634823 


1.992979 


2.428402 


2.003678 


2.634823 


2.970263 




20 


1.618424 


1.988671 


2.405376 


2.012146 


2.618424 


2.956329 


30 


50 


1.618036 


1.988497 


2.404869 


2.013483 


2.618036 


2.956009 




75 


1.618034 


1.988489 


2.404867 


2.013604 


2.618034 


2.956007 




100 


1.618034 


1.988486 


2.404867 


2.013644 


2.618034 


2.956007 


oo 


oo 


1.618034 


1.989044 


2.404867 


2.000000 


2.618034 


2.956295 



(mass 2mi), respectively), in agreement with results obtained by a truncated 
fermionic space approach. Clearly, our results also agree with the previous 
numerical investigations of the Ising model in a magnetic field, 8-12 ^}iich 
however are much less accurate. One disagreement is due to an error: Given 
the order of masses in equation (1.10) of Sagdeev and Zamolodchikov, ® the 
critical limit (their equation (1.9)) should read Aq = 0, Ai = 1/16, A2 = 1/2, 
A3 = 17/16, A4 = 2. 



7 Concluding Remarks 

The Bethe ansatz solutions of the off-critical dilute A3 model were studied 
numerically. Our results support the equivalence to Zamolodchikov's magnetic 
Ising model ^'■^ as well as recent assumptions on the string type of the solutions 
which enter a thermodynamic Bethe ansatz calculation ^® of the ^'-matrix. A 
more detailed account of our results will be published elsewhere. ■^■^ 
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